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§ ■ ABSTRACT 

O , The variability properties of accreting compact objects predicted by dynamical 

| models are characterized by a number of distinct frequencies that are specific to 

each model. Because the accretion disks around neutron stars and black holes are 
hydrodynamic flows, these modulation frequencies cannot be strictly equal to test 
particle frequencies. I discuss several implications of hydrodynamic corrections to the 
modulation frequencies predicted by dynamical models. Finally, I show that the recent 
detection of yet a third kHz QPO in three neutron-star systems favors a previously 
developed model that attributes the various QPO frequencies to fundamental general 
relativistic frequencies in the accretion flow. 
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1. INTRODUCTION 

The accretion flows around compact objects are turbulent and are, therefore, expected to 
be highly variable. The frequencies that correspond to the dynamical timescales can be as high 
as ~ 1 kHz and depend strongly on the distance from the central object. For example, for a 
Kerr spacetime and in the limit of slow rotation, the Keplerian orbital frequency at a radius r 
scales as ~ r -3/2^ periastron precession frequency as ~ r - 5 / 2 5 whereas the nodal precession 
frequency (if it is caused by frame dragging) scales as ~ r~ 3 (see, e.g., Perez et al. 1997; Stella, 
Vietri, & Morsink 1999). Because of such strong dependences, it is surprising that the observed 
power-density spectra of accreting systems are characterized by a distinct number of frequencies 
in the form of relatively coherent {/S.v/u ;> 10 -2 ) quasi-periodic oscillation (QPO) peaks. 

Modeling theoretically the presence of multiple QPO peaks is challenging. Current models 
seek to identify characteristic radii in the accretion flow that pick only a small range of frequencies 
and hence produce narrow QPO peaks in the power spectra. Such radii include the radius of the 
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innermost stable circular orbit (which is important in the diskoseismic models of, e.g., Nowak & 
Wagoner 1991; Perez et al. 1997; Silbergleit & Wagoner 2000), the magnetospheric radius (Alpar 
& Shaham 1985, Strohmayer et al. 1996), or the sonic-point radius (Miller, Lamb, & Psaltis 1998). 
The simultaneous presence of low- and high-frequency QPOs is then attributed either to multiple 
characteristic radii in the accretion flow (as discussed, e.g., in Miller et al. 1998; Psaltis et al. 
1999b) or to different characteristic frequencies occurring at the same region in the flow, as in 
the case of the diskoseismic models (see Wagoner 1999) or the relativistic models of Stella et al. 
(1999) and Psaltis & Norman (2000). In the first class of models, the observed tight correlations 
between the different QPO frequencies (see, e.g., Psaltis et al. 1999b; Psaltis, Belloni, & van der 
Klis 1999a) may be accounted for by combining the dependence of each QPO frequency on the 
mass accretion rate (see, e.g., Psaltis et al. 1999b). In the latter case, frequency correlations can be 
easily reproduced since all QPO frequencies are determined almost entirely by a single parameter 
(e.g., the characteristic transition radius; Psaltis & Norman 2000). 

A characteristic signature of each QPO model is the spectrum of QPO frequencies that it 
predicts beyond the ones currently detected (see, e.g., Miller 2000). For this reason, the inadequacy 
of current data to distinguish between different models can be overcome with the detection of 
additional QPO peaks. Even though some attempts for detecting expected QPO peaks at different 
frequencies have given negative results (Mendez & van der Klis 2000), a third kHz QPO has been 
recently discovered in three atoll sources (Jonker, Mendez, Sz van der Klis 2000). This new QPO 
is significantly weaker than the other kHz QPOs and has a frequency that is ~ 60 Hz higher than 
that of the previously known lower kHz QPO. 

Motivated by this detection, I discuss in this Letter the importance of the hydrodynamic 
corrections to the predicted QPO frequencies that are introduced by the previously developed 
dynamical model that identifies the various QPO frequencies with general relativistic frequencies 
in the accretion disk (Psaltis & Norman 2000). I then use this property to show that the model 
can account for the observed properties of the newly discovered QPO. 



2. HYDRODYNAMIC EFFECTS ON THE PREDICTED QPO FREQUENCIES 

The high frequencies that characterize the observed QPOs in neutron star sources imply that 
the accretion flows around them must be modulated at the dynamical timescale. A simple estimate 
of the three characteristic dynamical timescales at distance r away from a compact object is given 
by the three test-particle frequencies, i.e., the orbital frequency in the azimuthal direction, the 
epicyclic frequency in the radial direction, and the oscillation frequency in the vertical direction. 
For a slowly rotating compact object these frequencies are (see, e.g., Perez et al. 1997) 

H 2 = — ^, (1) 

r 3 [l + a*(M/r) 3 / 2 ] ' V ; 



.2 _ n 2 



n 2 (l - 6M/r) , (2) 
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and 

n 2 ± = n 2 {l-4a*(M/rf/ 2 } , (3) 

respectively, where M is the mass and a* the specific angular momentum per unit mass of the 
compact object, and I have set the fundamental constants to c = G = 1. The resulting orbital 
frequency (Q/2tt), periastron precession frequency [($7 — k)/2it], and nodal precession frequency 
[(0, — Q±)/2ir] were shown to be comparable to and follow similar correlations with the three 
observed QPO frequencies (Stella & Vietri 1998, 1999; Stella et al. 1999). 

An accretion disk, however, is a hydrodynamic flow and hence the test-particle frequencies (|])- 
@ cannot be strictly equal to the frequencies at which the flow is modulated. For example, 
the frequencies of the trapped g— and c— modes in the model discussed in Wagoner (1999) 
are determined mostly by the maximum of the radial epicyclic frequency (eq. ||]) and the 
nodal-precession frequency (ft — fij_) but also depend weakly on the local sound speed. In the case 
of the model discussed in Psaltis & Norman (2000), the inner accretion disk is modulated mainly 
at three characteristic frequencies, 



2tt/i = n - n 



nl(nl-K 2 ) 



(4) 



2n 2 d n 2 ± + 2(k 2 - n\) 2 

27Tf 2 = n-K[ 1 + ^|) , (0) 



and 

2vr/ 3 = fl , (6) 

where O c = c s /r, c s is the speed of sound, and 0,^ ~ u r /r is the inverse radial drift timescale at the 
transition radius r. In the notation of Psaltis & Norman (2000), the mode (m = 1, n = 0) provides 
the biggest contribution to the peak at frequencies f 2 and /3, whereas the mode (m = l,n = 1) 
provides the biggest contribution to the peak at frequency f\. 

It is important to note here that the corrections to the three dynamical frequencies introduced 
by the model are specific to the fact that they describe the modulation frequencies of the density 
in the inner disk and depend on the details of the calculation. For example, because of the 
approximations used, the effects of radial pressure forces, which will affect the hydrodynamic 
corrections, have not been taken explicitly into account. The corrections introduced to the 
modulation of the X-ray luminosity will be different in detail but similar in magnitude, if such 
effects are taken into account. However, hydrodynamic corrections are generic to a dynamical 
model and have three important implications. 

First, the observed QPO frequencies will not be strictly equal to the general relativistic 
frequencies; the corrections can be as large as 10 — 30%, depending on the local sound speed 
and the viscosity in the accretion flow (see, e.g., Psaltis & Norman 2000). This property resolves 
entirely the small but statistically significant discrepancy in the fits of the uncorrected general 
relativistic frequencies to observations (Psaltis et al. 1999b; Markovic &; Lamb 2000). An example 
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Fig. 1. — The effect of hydrodynamic corrections on the frequency difference 73-/2, calculated 
for a non-rotating compact object of mass 2.05M Q , a viscosity parameter a = 0.2, a fractional 
width of the transition region of 5r/r = 0.01, and different values of the disk scale-height h/r (see 
Psaltis & Norman 2000 for definitions). The data points are for the neutron-star source Sco X-l 
and correspond to the peak separation of the kHz QPOs and the frequency of the upper QPO peak 
(after van der Klis et al. 1997). 

is shown in Figure 1, where the predicted frequency difference fz — ji is compared to the changing 
peak separation of the kHz QPOs in the neutron-star source Sco X-l (van der Klis et al. 1997), for 
standard values of the model parameters. The observed decrease of the kHz QPO peak separation 
with increasing frequency is slower than predicted by the model but can be accounted for with 
a simultaneous small increase of the disk scale height. The resulting frequencies are in good 
agreement with the data without the need for introducing characteristic frequencies of elliptical 
orbits (as originally suggested by Stella & Vietri 1998), which are unlikely to be sustained in an 
accretion disk. Note here that the results plotted in Figure 1 are for a non-rotating compact object 
and were not obtained from formally fitting the data, but are shown to illustrate the effect of 
hydrodynamic corrections. However, detailed comparisons with data of the dynamical frequencies 
calculated for realistic neutron-star spacetimes were shown to affect only slightly the conclusion 
(Markovic & Lamb 2000). 
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Second, the centroid frequencies of QPO peaks produced by resonances that correspond to 
harmonics of general relativistic frequencies are not exactly harmonically related. This is true, 
because the hydrodynamic corrections are different for the different modes and hence the ratio of 
the frequencies of the modes is not an integer number but depends on the properties of the flow. 
Note, however, that if the presence of additional harmonics besides the fundamental is caused 
by the non-sinusoidal profile of the modulation of the X-ray flux and not by the presence of 
higher-order modes, then the frequency ratio will be an integer. This might lead to the presence 
of multiple peaks at nearby frequencies with possibly different FWHMs, as for example observed 
in the case of Cyg X-l by Nowak (2000). 

Finally, different modes may contribute power to the same QPO peak but at slightly different 
centroid frequencies, again because of the different hydrodynamic corrections. For example, all 
the (m = l,n > 0) modes in the analysis of Psaltis & Norman (2000) contribute to the variability 
power at frequencies comparable to ~ /2, i.e., to the lower kHz QPO (the contribution of the 
modes with n > is small and was not explicitly pointed out in the original analysis, even though 
it was present). Indeed, expanding the general expression for the response function (eq. [27] in 
Psaltis & Norman 2000) for the (m = l,n = 1) mode and at frequencies 2tt/ ~ fl — k, it simplifies 
to 
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when (5r /r)(K/fl c ) 2 <C 1. Equation (0) describes a narrow resonant mode centered at a frequency 
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where = ( r /Sr)fl^ and 5r/r is the fractional width of the transition region. This resonance is 
to be compared to the corresponding resonance of the (m = l,n = 0) mode, which is described by 
(see, eq. [35] of Psaltis & Norman 2000) 
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when (5r/r) {n/fl rmc Y < l. 
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Fig. 2. — The response function for modes (to = 1, n = 0) and (to = l,n = 1) at high frequencies, 
calculated for 0/2tt = 1138 kHz, k/2tt = 280 Hz, (O - n ± )/2ir = 50 Hz, Q. c /2-k = 200 Hz > and 
L0\,/2Tr = 50 Hz. The parameters correspond to typical values of the hydrodynamic corrections in 
an accretion disk (see eq. [14]) and were chosen to reproduce the observed power spectrum of the 
source 4U 1728-34 (Jonker et al. 2000). 

Figure 2 shows the resulting response of these lowest-order modes at high frequencies. As it 
is also apparent from equations (||) and the two peaks at frequencies comparable to ~ /2 are 
displaced relative to each other and, in the particular model considered here, their separation is 
equal to 



As a result, two QPOs at frequencies comparable to that of the lower kHz QPO may appear in 
the power spectra of accreting compact objects. Note here that the two frequencies are identical 
in the limit £l c — > 0, i.e., for a purely kinematic model. 

The amplitudes of the two QPOs will depend on the power-spectrum of the driving 
perturbations, which might be different in the two modes that produce the two QPOs, even 
though their frequencies are very similar. (The maximum of the response function does not 
provide an estimate of the amplitudes of the QPOs.) However, even for a comparable amplitude 
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of perturbations in both modes, the resonance described by Aiq will produce the dominant QPO 
peak (see also Psaltis &: Norman 2000), whereas the resonance described by An will produce a 
weaker sideband, as Figure 2 suggests. In the absence of any additional broadening mechanisms, 
the FWHM of the two lower kHz QPOs (at frequencies f 2 and /2b) will be significantly smaller 
than the FWHM of the upper kHz QPO (at frequency /s). However, the two lower kHz QPOs 
may have a difference of up to ~ {Q±/k) 2 ^ 15 in their relative FWHM. 

3. DISCUSSION 

In the previous section, I discussed a number of implications of the hydrodynamic corrections 
introduced to the predicted QPO frequencies by dynamical models. Moreover, I showed that the 
previously developed model that attributes the observed QPOs to general relativistic frequencies 
in the accretion flow can account for the presence of multiple kHz QPOs, as recently observed 
by Jonker et al. (2000). In this section I will concentrate to the case of the well studied source 
4U 1728—34 (Jonker et al. 2000) and discuss the implications of the newly observed QPOs for 
various models. In this source, four QPOs are observed simultaneously at the following frequencies: 
z/i = 41.5 ± 0.2 Hz (the low-frequency QPO), u 2 = 795 Hz (the lower kHz QPO), u 2b = 859 ± 2 Hz 
(the newly discovered QPO), and ^3 = 1138 ± 2 Hz (the upper kHz QPO). 

According to beat-frequency models (see, e.g., Strohmayer et al. 1996; Miller et al. 1998), 
the highest QPO frequency ^3 is that of the Keplerian frequency of a characteristic radius in 
the accretion disk and the difference z/3 — v 2 is comparable to the spin frequency of the neutron 
star. In fact, oscillations during Type I X-ray bursts have been detected from 4U 1728—34 at a 
frequency of ~ 363 Hz, which are also thought to correspond to the spin frequency of the neutron 
star (Strohmayer et al. 1996). In this picture, the presence of the newly discovered QPO at a 
frequency that is ~ 64 Hz higher than v 2 might be accounted for if the radiation that emerges 
from the neutron star is modulated not only at the stellar spin frequency, but also at a frequency 
equal to ~ 64 Hz, introducing sidebands to the stellar spin (see, however, Alpar 1986 for an 
alternative way of generating sidebands in a beat-frequency model, because of the radial drift of 
the accreting material). As pointed out by Jonker et al. (2000), this latter frequency is not equal 
to the frequency of the observed QPO at v\ and therefore the accretion flow needs to produce 
modulations at two distinct, unrelated, low frequencies. One such low frequency is the nodal 
precession frequency (f2 — £l±)/2ir. Note, however, that, as also discussed by Psaltis et al. (1999b) 
for other similar cases, the QPO at frequency v\ cannot be caused by frame-dragging effects if the 
spin-frequency of the neutron star is equal to ~ 363 Hz. (In the relativistic model described in §2 
the frequency v\ is related to the frame-dragging frequency, but the peak separation of the kHz 
QPOs is not equal to the spin frequency of the neutron star.) This can also be extended to the 
case of the ~ 64 Hz frequency, which is even higher than v\. Indeed, such a high nodal precession 
frequency would require an unrealistically high neutron-star moment of inertia of (see, e.g., Psaltis 
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et al. 1999b) 



'* ^(^z)(M~ f^] (13) 



M/M & ^ ' \64 Hz/ V363 Hz/ V H38 Hz, 
where I45 = I /10 45 gem 3 is the moment of inertia, M/Mq the mass in solar units, and f sp i n the 
spin frequency of the compact object. Even if a frequency of ~ 64 Hz corresponded to the first 
overtone of the nodal precession frequency, the required moment of inertia / (M / Mq) ;> 1.55] 
would be higher than allowed by general relativity for a slowly rotating star, if any of the current 
equations of state for neutron-star matter were valid up to ~ 1.5 times the nuclear saturation 
density. As a result, none of the two required low frequencies can be equal to the frame-dragging 
frequency, if the spin frequency of the neutron star in 4U 1728—34 is equal to the burst oscillation 
frequency. 

In the case of the dynamical model described in Psaltis &; Norman (2000), the presence 
of multiple QPO peaks in the power-density spectra is a direct consequence of the fact that in 
strong-field general relativity the three dynamical frequencies (eq. [||-||3|) are significantly different 
from each other. In particular, all the QPO peaks observed so far can be accounted for by only 
the two lowest-order modes (to = 1, n = 0) and (to = 1, n = 1). 

In this interpretation, the difference /2b — H depends mostly on the radial epicyclic frequency 
and the inverse sound-crossing time (see, eq. [Q). For a standard a-disk, the latter is equal to 
£l c ~ (h/r)£l±, where (h/r) is the scale-height of the disk (see, e.g., Psaltis <fc Norman 2000). As a 
result, the predicted frequency separation between the two lower-kHz QPO peaks becomes 

{ ^ h ^(^(^\ 2 (^ r y^ (14) 



,0.127 V 1138 Hz/ V 280 Hz 
and hence the observed frequencies can be easily accounted for, for typical values of the model 
parameters (see also Fig. 2). 

Higher-order modes will also contribute to the variability power at frequencies comparable to 
/2 and may produce additional structure to the lower kHz QPO peak. The overall spectrum of 
observed modes will depend mostly on the amplitude of the perturbations that drive the variability 
outside the characteristic transition radius. A detailed analysis of the various predicted resonant 
frequencies and their amplitudes is beyond the scope of this paper and will be presented elsewhere. 
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